Ancilla-assisted sequential approximation of nonlocal unitary operations 
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We consider the recently proposed "no-go" theorem of Lamata et al [Phys. Rev. Lett. 101, 180506 
(2008)] on the impossibility of sequential implementation of global unitary operations with the aid 
of an itinerant ancillary system and view the claim within the language of Kraus representation. 
By virtue of an extremely useful tool for analyzing entanglement properties of quantum operations, 
namely, operator-Schmidt decomposition, we provide alternative proof to the no-go theorem and 
also study the role of initial correlations between the qubits and ancilla in sequential preparation 
of unitary entanglers. Despite the negative response from the no-go theorem, we demonstrate ex- 
plicitly how the matrix-product operator (MPO) formalism provides a flexible structure to develop 
protocols for sequential implementation of such entanglers with an optimal fidelity. The proposed 
numerical technique, which we call variational matrix-product operator (VMPO), offers a compu- 
tationally efficient tool for characterizing the "globalness" and entangling capabilities of nonlocal 
unitary operations. 
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I. INTRODUCTION 

Engineering arbitrary global unitary operations entan- 
gling simultaneously multiple qubits is generically re- 
garded to be a task of formidable difficulty, as it may 
require a complicated combination of exponentially many 
gates to approximate [l|. However, efficient implementa- 
tion of such genuinely entangling multiqubit operations is 
of paramount importance in quantum computation. Re- 
cently, it has been wondered if the situation could be fa- 
cilitated by devising an ancilla-assisted decomposition of 
a multiqubit unitary operation into a sequence of ancilla- 
qubit unitary operations where each qubit is allowed to 
interact locally and only once with an itinerant ancil- 
lary system (e.g., a trapped multilevel atom coupled to 
a single mode of an optical cavity in the realm of cavity 
QED experiments [2-4]) and without measurements [5| 
(see Fig. [T]) . Equivalently, according to the suggested sce- 
nario, instead of entangling directly neighboring qubits k 
and k + 1, an ancilla a intervenes to get entangled with 
qubit k and the ancilla state is swapped afterward with 
that of qubit k + 1. Such decomposition, however, im- 
mediately faced a resounding no and was proved to be 
impossible in general for genuinely entangling unitaries, 
with the controlled-NOT (CNOT) gate being the hand- 
iest counterexample [5j. A reductio ad absurdum proof 
strategy for such a "no-go" theorem implied the incom- 
patibility of unitarity of the last constituent two-qubit 
operation and the deterministic nature of the protocol. 
More precisely, having introduced an ancillary system to 
convey entanglement throughout the register of qubits, 
we eventually wish the ancilla decouples from the qubits 
in the very last step. But it turns out that for arbitrary 
entangling unitary the ancilla always remains entangled 



with the qubit system, thereby spoiling the promised de- 
terminism. 

In the present paper, we reconsider the issue of such a 
sequential quantum "factory" of unitary entanglers within 
the framework of Kraus representation Q and operator- 
Schmidt decomposition 0, H| and provide alternative 
proofs to the no-go theorem forbidding the realization of 
such a setup. It turns out that employing such tools pro- 
vides a systematic and transparent approach to various 
aspects of the problem and reveals the insurmountable 
mathematical obstacles to realization of such a scenario. 
Moreover, despite this negative answer, we investigate 
the possibility of implementing a nonlocal unitary with 



(a) 



(b) 



U 



12- Na 



u 



Na 



u 2a u la 



a 
N 

2 
1 







i 


( 


s 



a 
N 

2 
1 




-o- 



Figure 1. Quantum circuit representation of (a) a generic non- 
local multiqubit unitary and (b) its ancilla-assisted sequential 
implementation. Each bipartite unitary Uka acts only on the 
Hilbert space of qubit k and ancilla a and leaves other qubits 
unchanged. To emphasize this fact, instead of the common 
circuit representation of an arbitrary bipartite unitary by a 
rectangular box, we have used vertical lines with open circles 
at each qubit they act upon. Ancilla states are shown with 
solid bold lines throughout. 
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a certain imperfect fidelity. In other words, we demon- 
strate it is always possible to satisfy both conditions of 
the sequentiality and unitarity of the ancilla-qubit opera- 
tions at the same time, however, at the price of ending up 
with a sequentially implemented unitary whose action is 
not perfectly equivalent to the original global unitary but 
is rather closest to that in some sense. This will be real- 
ized by exploiting the tools from matrix-product operator 
(MPO) theory \mT\\. We also study the role of initial 
correlations between the ancilla and qubits upon sequen- 
tial preparation of unitaries and the way they could ever 
affect the evolution of the joint ancilla-qubit system. 



II. SEQUENTIAL QUANTUM "FACTORY" OF 
UNITARY ENTANGLERS VIEWED IN THE 
LANGUAGE OF QUANTUM OPERATIONS 

Let us first illustrate explicitly what is exactly meant 
by claiming that no entangling unitary can be imple- 
mented in a sequential way: For this purpose, we suggest 
viewing the problem within the framework of quantum 
operations [TJ, |6( where the qubit chain is regarded as the 
principal system and the ancilla as the environment (see 
Fig- HI- For simplicity, let us consider for now the specific 
case of only two qubits, N — 2. We shall later general- 
ize our argument to the case of an arbitrary number of 
qubits N. Assuming that the qubits and ancilla start out 
in a product state /?i 2a — P12 ® Pa, the evolution of the 
state of qubits is given by 

£{pu) = Tr a [C/ 12a (p 12 ® p a )Ul 2a ] = E k pi2E\ , (1) 

k 

where Ek = a (k\Ui2 a \0) a are Kraus operators in an 
operator-sum representation [l[ of quantum operation £ 
and \k) a represent an orthonormal basis for the (finite- 
dimensional) state space of the ancilla; in particular, |0) 
is assumed to be the initial state of the ancilla. We note 
that the aim is here to implement an entangling target 
two-qubit unitary U\2 with the aid of an ancilla a whose 
action should be eventually factorized as U\2a = £/i2® la- 
Should the latter be the case, the Kraus operators are 
simply given by 



E k . 



|tfl2®la|0>„ = U 12 S k0 



(2) 



The sequential preparation of the two-qubit unitary U12, 
on the other hand, suggests Kraus operators of the form 



E k = a (k\{U 2a ® li)(l 2 ® U la )\0) a 



(3) 



Note that at step k of the sequential implementation of 
a global A^-body unitary operation U\2...n, the two-body 
unitary Uka entangles only the ancilla a and qubit k and 
leaves other qubits unchanged, the latter action having 
been denoted by a properly indexed identity operator 1 
above. Now, if the entangling unitary [7 12 could be im- 
plemented sequentially, E\ in Eq. ^ should also vanish. 
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Figure 2. Ancilla-assisted implementation of entanglers 
viewed as a quantum operation £(pi2...n)- The joint ancilla- 
qubit system evolves under a unitary transformation Ui2...Na- 



Acting Ei on some initial state of qubits 12 and per- 
forming Schmidt decomposition on the partition l\a of 
the entangled state Ui a \00) la yield 

£l|0z) 12 = a(l\{U 2 a ® ll)(l 2 <8> E/la)|(M)}i 2a 

X 

= ^ a (l|C/2aNj) 2 a|^} 1 =0, (4) 
3=0 

with the Schmidt rank \ > 1 fl2l |. and by convention 
we have incorporated the Schmidt coefficients into the 
Schmidt vectors [l3[. The orthogonality (and linear in- 
dependence) of \4>j) 1 implies that the only way for E\ to 
vanish for a sequentially prepared unitary would be that 
o(l|C^2a|ii)2a vanish for all i and j = 0, 1, ... . The lat- 
ter, however, causes two rows of £/ 2o and, consequently, 
its determinant to vanish, thereby violating the invert- 
ibility and, in turn, unitarity of Z7 2q . 

For sequential implementation of a global unitary one 
should satisfy two constraints of sequentiality and unitar- 
ity of the constituent ancilla-qubit operations. However, 
as illustrated above, these two cannot be fulfilled at the 
same time. 

Note that the result above may be readily extended 
to an arbitrary number of qubits A^ since it is always 
possible to map this general case to that of only two 
unitaries U^aUia by combining the first n — 1 unitaries 
and the rtth one, where n is the smallest step at which the 
ancilla decouples from the qubit chain [5j- Therefore, for 
the sake of brevity, we shall henceforth only consider the 
paradigmatic case of A^ = 2 and keep in mind that the 
results will accordingly be valid for an arbitrary number 
of qubits A^ too. 

We point out the "no-go" statement above holds even if 
we relax the assumption that the ancilla and qubits start 
out in a product state. To see this, we exploit the Fano 
form of the density matrix (l4l . Qj| to follow the evolution 
of the joint ancilla-qubit system in the presence of initial 
correlations between the qubits as the open system (l6| 
and the ancilla as the environment. Assuming the action 
of the ancilla can be factorized like U\2 a — U12 ® 1 Q , the 
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quantum operation then is given by [17j 

£(Pl2) = U 12 pi2U{ 2 

+ J2«( k \ U ^l>l 2 ® T?Ill 2a \k) a , (5) 
fe 

where <jj and Tj are the generators of SU(4) and SU(D), 
respectively, and D is the dimension of the Hilbert space 
of the ancilla. Straightforward algebra then yields 

£( P i2) = U 12 p 12 U\ 2 + ]T U 12 al 2 ul 2l ^ a T r {Tf} . (6) 

id 

But generators of SU(-D) are known to be traceless. So 
we conclude that for the case where the evolution of the 
whole ancilla-qubits is factorable with respect to the an- 
cilla (or the ancilla decouples from the qubit chain), the 
initial correlations would not affect the final states of the 
qubits. Contrarily, for a nonfactorable sequentially im- 
plemented unitary, the correlation term in Eq. ^ does 
not necessarily vanish, and initial correlations may play a 
significant role in the evolution of the joint ancilla-qubit 
system and, in turn, in the final states of both qubits and 
the ancilla. 

The main obstacle in sequential implementation of 
an entangling unitary, namely, the fact that the ancilla 
cannot be set to decouple from the qubit chain in the 
last step, may also be expressed in terms of the well- 
known problem of separability criteria in quantum infor- 
mation [la - |2lft : When we start from a 1-qubit bisepa- 
rable state [22J of ancilla and qubits pi 2a = p\ 2 ® p a , 
straightforward algebra employing Kraus representation 
for both the ancilla and qubits implies that a factorable 
evolution U\ 2a = U\ 2 ® U a leaves such separability in- 
variant, i.e., p' 12a — p\ 2 ® p' a where a prime indicates the 
final states. However, this ceases to be the case for the 
ancilla-assisted sequential evolution considered above. 



III. ALTERNATIVE PROOF TO THE NO-GO 
THEOREM BASED ON OPERATOR-SCHMIDT 
DECOMPOSITION 



In this section we provide an alternative proof to 
the no-go theorem by virtue of an extremely useful 
and elegant decomposition known in operator algebra as 
operator- Schmidt decomposition [7, 8] according to which 
the bipartite ancilla-assisted unitaries may be decom- 
posed as 



U 2a = J^A 2 ®B?, 



(7a) 
(7b) 



with {A?}, {A}}, {Bf }, and {Cf} being some orthogo- 
nal (but not necessarily normalized) operator bases de- 
fined on the Hilbert-Schmidt space of operators acting on 



qubits 2 and 1 and ancilla a, respectively. Of particular 
importance is \2 (xi)i called the Schmidt number of op- 
erator U 2a (Uia), defined as the number of nonzero coeffi- 
cients in its operator-Schmidt decomposition. Applying 
the decomposition to each bipartite entangling unitary 
U ka yields 

U 12a = (U"aa®ll)(la®^la) 

= (J2 A] ® B\ ® li)(l 2 ® A ) ® C V 



Y,A l j ®A 2 ®{B«C?) 



(8) 



Without loss of generality we may take \2 — 2; the 
same argument would hold for an entangling unitary of 
Schmidt number \2 = 4 [2^|. In order that the ancilla 
decouples unitarily from the qubit chain, it should be 
possible for it to be factorized like BfC" = U a for all i 
and j and with U a some unitary operator on ancilla space. 
This in particular implies that BfCf = B 2 C^, eventually 
suggesting U 2a — (X^j A?) ® B a to be a nonentangling 
operation. The latter, however, is in contradiction to 
the original assumption that U 2a has a Schmidt number 
larger than 1, and this readily completes the alternative 
proof. 



IV. OPTIMAL IMPLEMENTATION OF 
NONLOCAL UNITARIES WITHIN THE 
FRAMEWORK OF MATRIX-PRODUCT 
OPERATORS 



Having accepted the impossibility of implementing a 
nonlocal unitary in a sequential way, we may wonder if 
we could rather try to approach a target unitary entangler 
within a sequential procedure, albeit in an approximate 
way and with an optimal fidelity. More precisely, the 
aim is then to find bipartite unitaries U\ a and U 2a that, 
when applied sequentially instead of the global unitary 
U\ 2a , minimize the Frobenius norm distance [24] 



C 



^13o-(^2o®ll)(la®^lo) 

2D - 2Re{Tr[t/ 1 t 2a (£/ 2a ® 1 1 )(1 2 ® U la )] } . (9) 



Note here that we hereby require both sequentiality and 
unitarity conditions to be fulfilled without violating the 
above-mentioned no-go theorem: In fact, these two plus 
the implicit condition that the sequentially implemented 
unitary should produce the same action of the global one 
(i.e., the identity U 12a \^) = {U 2a ® li)(l 2 <S> U la )\V) for 
arbitrary state |^)) do save the no-go theorem and make 
such a distance never vanish. Without such a condition, 
we could, of course, freely choose the sequential unitaries 
so that the ancilla decouples unitarily in the last step. 

The constrained nonlinear optimization problem (25| 
associated with the cost function above may be solved ef- 
ficiently by a variational search in the space of the MPOs 



4 



h 



U 



[i] 



h 



J2 



u 



[2] 



Jn 



u 



[N] 



l N 



u 



[12-JV] 



72 



Jn 



arising from the sequential nature of theprotocol. To see 
how operators of matrix-product form [9l-lll| arise natu- 
rally within such a prescription, we express each ancilla- 
qubit unitary in terms of an orthonormal basis 
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(10) 
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Figure 3. Graphical representation of the overlap of MPOs 
used to calculate the contractions in the cost function 
Eq. (|12[) . The square boxes represent U< ka i matrices of a 
sequentially prepared unitary, and the rectangular box rep- 
resents the target global unitary. The links connecting the 
boxes represent indices that are being contracted (or summed 
over) . 
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where Roman (Greek) letters denote qubit (ancilla) in- 
dices. The sequentially implemented unitary, as a result, 
reads 



12 ,jz 0.2 ,@2 h ,jl ai ,Pi i-2 ,]2 ii ,ji 



l % 32 U^\i 2 )(h\®\h)(h\ 



(11) 



r 



where we have defined in the last line 
Eatft ^jtlh \ a k){Pk I = U^f" to cast the original 
expression into products of unitary operators. We 
may now easily recognize (fTTj) as an MPO already 
encountered in the calculation of finite-temperature 
density matrices (26j. With this, the cost function (|9]) in 
the MPO representation reads 



C = 2D 

- 2Re { Tr (E E E u w 3 ^ u m nu w n 

f 2 J'l '2J2 HJl 

®|jaJi>(?2ji|)} . 



(12) 



where we have expressed the target unitary too, in terms 
of an orthonormal basis, and have accordingly defined 



E^ j2 'v></3'i 

ol' 



[21] 



to have 



^ = EE^f'^i^i)02iii 



(13) 



(14) 



The calculation of the cumbersome summations on the 
right-hand side of Eq. (fT2j) finds a very simple form 
in MPO representation as illustrated in Fig. [3] The 
minimization problem can be done very efficiently in 
an iterative variational optimization technique (similar 



to those already used in the context of matrix-product 
states [3, [27l - [30| ) in which we fix all but one of the se- 
quential unitaries, let's say U\u\, and minimize the cost 
function Eq. (|12[) by varying over the matrix elements of 
U\ k -\ . In the next iteration, the neighboring unitary is op- 
timized, and once all the unitaries so obtained have been 
optimized locally, we sweep back again and so on until 
convergence. In that sense, we suggest an extension of 
the previously developed method of variational matrix- 
product state (VMPS) (2l413| from states to operators 
and suggest calling it the variational matrix-product op- 
erator (VMPO) method. 

In the case that the ancilla is a qubit, i.e., D = 2, the 
two-qubit sequential unitaries of the variational MPO in 
Eq. (jlip may be expanded in terms of the complete basis 
of Pauli matrices: 



Ik] 



exp(- 



3 

; E 

Jl J2=0 



(15) 



where ,- 2 are real- valued coefficients and a±, 02, and (73 
are the usual Pauli sigma matrices, with tro = 1 being the 
identity matrix. The minimization of the cost function 
Eq. p2p then amounts to finding the optimal coupling 
matrix h\ k ^ at each step of the iterative optimization. The 
fidelity of the sequential implementation of the global 
unitary U\2 a may be quantified then as J- = 1 — C, where 
C denotes the normalized converged cost function C of 
Eq. (I12p upon the very last sweep. The normalization is 
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Table I. The values of the Frobenius fidelity for sequential 
implementation of various global gates with an ancilla of di- 
mension D — 2. 



Unitaries: CNOT CZ CPHASE SWAP Toffoli Fredkin 
Fidelity 0.7071 0.7071 09239 O50 075 075 - 



taken care of after dividing by the Frobenius norm of the 
involved operators. 

We have applied the outlined procedure to some 
paradigmatic gates of quantum computing when the an- 
cilla has the dimension D = 2, as illustrated in Ta- 
ble ID For instance, the two locally equivalent gates 
controlled-NOT (CNOT) and controlled-Z (CZ) can be 
implemented sequentially with a fidelity of 70.71%. The 
29.29% error we associate with the existence of a "fi- 
delity gap," a strict theoretical "red line" beyond which 
the global unitary cannot be implemented in a sequential 
fashion. 

Similar optimization techniques may be exploited for 
higher ancilla dimension D = 4, however, upon using the 
generators of SU(4) for the ancilla instead of Pauli matri- 
ces in Eq. (|15[) to see if it has any influence on the values 
of the gaps. We have performed such a simulation for the 
same gates as in Tableland have found that the values of 
the gaps remain unchanged for D = 4 compared to that 
of D = 2. Since D — 4 is the maximum possible ancilla 
dimension for both N — 2 and N — 3 cases [3l| , we con- 
clude that the reported values of the gaps are intrinsic 
to the gates, irrespective of the ancilla dimension. 

We believe also that such intrinsic values of the gaps 
are closely related to the entanglement content of each 
unitary operation, the figure of merit for this purpose 
being the Schmidt strength [8] (proportional to the sum 
of nonzero Schmidt numbers of an operator) as a mea- 
sure of operator entanglement [32|-[34[ ■ Exemplifying val- 
ues of Schmidt strengths are 0.5 and 0.75 for CNOT and 
SWAP, respectively. As expected, the larger the Schmidt 
strength of the gate is, the larger the corresponding fi- 
delity gap is. 



V. CONCLUSIONS 

In conclusion, we have studied the problem of se- 
quential implementation of entangling unitary operations 
within the various languages of Kraus representation, 
operator-Schmidt decomposition and MPO representa- 
tion. Exploiting such tools allowed us to revisit the previ- 
ously proposed no-go theorem through a systematic and 
instructive scheme that reveals various mathematical as- 



pects of the subject. The Kraus representation approach, 
in particular, enabled us to probe possible correlations 
between qubits and an ancilla after unitary evolution of 
the joint system and thereby to compare the ultimate ac- 
tion of sequentially prepared unitary to that of a global 
entangler. Employing MPO formalism, we also devel- 
oped numerical techniques for an efficient realization of 
a sequential version of an entangling unitary operation 
with an optimal fidelity. As a realistic scenario, the recipe 
could instruct an experimentalist how to obtain optimal 
two-qubit operations that lead to a sequential version of 
the desired multiqubit entangler with the highest possi- 
ble fidelity and without the need to perform any mea- 
surement. 

We stress that these results and the proposed optimiza- 
tion protocols can be of wide potential applicability at the 
interface of condensed matter physics, quantum optics, 
and quantum information. In particular, they will be of 
direct relevance to any sequential physical setup like pho- 
tonic qubits, superconducting qubits, or quantum dots. 
From a methodological point of view, the proposed nu- 
merical technique, the VMPO method, promises to be a 
computationally efficient tool for variational search into 
the Hilbert-Schmidt space of operators. We believe that 
the method may particularly be utilized for characteriz- 
ing the entangling (and disentangling) properties of mul- 
tiqubit unitary operations. 

We point out, finally, that the no-go theorem does not 
refrain from implementing an entangling unitary if (in- 
stead of targeting the nonlocal unitary within a single 
sweep of the itinerant ancilla along the qubits chain) we 
allow the ancilla to sweep back and forth through the 
qubit chain and thereby perform multiple bipartite op- 
erations at various stages on each qubit. For instance, 
it can be easily verified that CNOT can be implemented 
through two rounds of ancilla operations each performed 
in a sequential manner. 
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